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Abstract 

The constrained Next-to-Minimal Supersymmetric Standard Model (cNMSSM) 
with mSugra-like boundary conditions at the GUT scale implies a singlino-like LSP 
with a mass just a few GeV below a stau NLSP. Hence, most of the squark/gluino 
decay cascades contain two r leptons. The gluino mass > 1.2 TeV is somewhat larger 
than the squark masses of > 1 TeV. We simulate signal and background events for 
such a scenario at the LHC, and propose cuts on the transverse momenta of two 
jets, the missing transverse energy and the transverse momentum of a hadronically 
decaying r lepton. This dedicated analysis allows to improve on the results of generic 
supersymmetry searches for a large part of the parameter space of the cNMSSM. The 
distribution of the effective mass and the signal rate provide sensitivity to distinguish 
the cNMSSM from the constrained Minimal Supersymmetric Standard Model in the 
stau-coannihilation region. 



1 Introduction 



The Next-to-Minimal Supersymmetric Standard Model (NMSSM, for recent reviews see [TJ 
Ej) provides an elegant solution to the /i-problem of the Minimal Supersymmetric Standard 
Model (MSSM) [3]: after the replacement of the /x-term in the superpotential of the MSSM 
by the coupling to a gauge singlet superfield S, the superpotential is scale invariant and the 
only dimensionful parameters in the Lagrangian are soft Supersymmetry (SUSY) breaking 
terms. 

After electroweak symmetry breaking, all components of the gauge singlet superfield S 
mix with the components of the MSSM-like Higgs doublet superfields H u and Hd- Accord- 
ingly NMSSM specific phenomena can take place in the CP-even Higgs sector, the CP-odd 
Higgs sector and the neutralino sector. In this paper we concentrate on the neutralino 
sector, where the 5th singlet-like neutralino (in addition to the two neutral higgsinos, the 
neutral wino and the bino) can be the Lightest Supersymmetric Particle (LSP) [1HS] and, si- 
multanously, give rise to a dark matter density in agreement with WM AP constraints [SHE] ■ 

Such a scenario is not far fetched: A conceptually simple origin of soft SUSY breaking 
terms is a minimal coupling to supergravity (with a flavour independent Kahler potential 
and minimal gauge kinetic terms), and the assumption of spontaneously broken SUSY in 
a hidden sector. Then, the soft SUSY breaking terms are universal at the Planck scale 
(not far from the GUT scale) in the form of universal scalar masses m , universal gaugino 
masses M\n and universal trilinear scalar couplings proportional to Aq. 

The correspondingly constrained NMSSM (cNMSSM) has been studied first in [T5] - [20] . 
Since then, the precision of the radiative corrections has been considerably improved, and 
the computation of the dark matter relic density has become possible [9J. Imposing the 
requirement of a dark matter relic density in agreement with WMAP constraints [2] as 
well as present constraints on Higgs and supersymmetric particle (sparticle) masses, the 
parameter space of the cNMSSM has been analysed recently in [2T1[22]. (For studies of 
the semi-constrained NMSSM, where the singlet- dependent SUSY breaking parameters are 
allowed to be non- universal, see [r2 ^ ir5 }l2"5H2"5] .) 

Within the cNMSSM, the dark matter constraints require a singlino-like LSP. The origin 
of this phenomenon is quite easy to understand: first, the CP-even scalar singlet s has to 
assume a non-vanishing vacuum expectation value (vev) in order to generate the required 
/x-term. For this reason its SUSY breaking mass squared m| at the electroweak scale must 
not be very large; otherwise the minimum of its potential is at s = 0. Second, m| is hardly 
renormalized between the GUT and the electroweak scales, which leads to m| ~ with 
the consequence that, in the cNMSSM, mo must be small compared to M\/2 and A$. 

It is well known that, within the cMSSM [26], a small value of mo would imply a stau 
(ti) LSP, which is not a reasonable candidate for the dark matter. In the cNMSSM, the 
singlino-like neutralino Xi can ~~ an d must - be lighter than the r± for this reason. In order 
for its dark matter relic density not being too large, its mass must be only a few GeV below 
the mass of the T\ such that x\ —T\ coannihilation processes are sufficiently fast [T2" | |2"T } |2"2"] . 
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Such a scenario would necessarily have a strong impact on sparticle searches at colliders 
(see [271125] for reviews of searches at the LHC): since the T\ is the NLSP and the singlino- 
like neutralino Xi couples only very weakly to all sparticles, the sparticles decay first into 
the ri under the emission of at least one r-lepton. Subsequently the t\ decays into Xi + r ) 
hence each sparticle decay cascade contains typically two r-leptons. In the present paper 
we study, for the first time, the corresponding implications for sparticle searches at the 
LHC. 

As discussed in [2~T||2"2"]. the sparticle (and Higgs) spectrum is quite constrained in the 
cNMSSM, and can essentially be parametrized by M\ii- The present lower bounds on 
m^ x and the lower LEP bound of ~ 114 GeV on the CP-even Standard Model like Higgs 
mass require Myz > 520 GeV which, inspite of m < 50 GeV, implies a quite heavy 
sparticle spectrum: squark masses > 1 TeV (apart from the somewhat lighter stop and 
sbottom masses), and a gluino mass > 1.2 TeV. Since consequently squark production will 
be relatively frequent compared to gluino production (in fact, squark + gluino production 
dominates), the signal events will mostly lead to two jets (not counting ISR and FSR) with 
a quite large as well as a large missing transverse energy E™ ms . The sum of the squark 
and/or gluino production cross sections is, however, just ~ 1 pb or less (for larger Mi/ 2 ). 

Each sparticle production event in the cNMSSM will contain typically four r-leptons in 
the final state: two of the r-leptons (those originating from t\ decays into Xi + r ) wm be 
quite soft due to the small t\ — Xi mass difference < 5 GeV and are difficult to detect, 
whereas the r-leptons from decays as x\ ~^ + t are relatively energetic. 

The aim of the present paper is to show that a dedicated analysis allows for much 
better signal to background ratios for the cNMSSM than standard (generic) supersymmetric 
analyses. We simulate and study signals and various Standard Model backgrounds for the 
LHC at 14 TeV cm. energy, and find that the signal to background ratio is sufficiently 
large allowing for the discovery of the cNMSSM for a wide range of values of M 1 / 2 . 

Events with E™ 1SS , jets and r-leptons could also be a signal for the (c)MSSM in the 
so-called stau-coannihilation region [2"5H3"3"] . Since the squark/neutralino spectrum of the 
cMSSM is necessarily different from the cNMSSM, the combination of M e ff (essentially the 
sum of the transverse momenta and E™ 1SS ) and the signal rate has sensitivity to distinguish 
the two models. 

In the next section we present the spectrum of the cNMSSM. In section 3 we discuss the 
signal, backgrounds and appropriate cuts. Details of the simulation of the cNMSSM signal 
and the Standard Model backgrounds, the effect of cuts, and the resulting cross sections and 
signal to background ratios are given in section 4. Section 5 is devoted to the comparison 
of the cNMSSM with the cMSSM, and Section 6 to conclusions and an outlook. 

2 The spectrum of the cNMSSM 

The NMSSM with a scale invariant superpotential W [H[2j differs from the MSSM through 
the replacement of the \i term in Wmssm by the coupling to a gauge singlet superfield S 
and a trilinear S self-coupling: 

Wmssm = ^H u H d +...-> W NM ssm = \SH u H d + ^S 3 + . . . , (2.1) 
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where we have omitted the quark/lepton Yukawa couplings. Hence, if the vev s of S is 
non-zero (induced by the soft SUSY breaking terms), an effective \i term /i e yj = As of the 
desired order of magnitude is generated and solves the \x problem of the MSSM [3J. 

Apart from the Yukawa couplings in the superpotential and the gauge couplings, the 
Lagrangian of the NMSSM depends on soft SUSY breaking scalar masses for the Higgs 
fields H u , Hj and S, the squarks and the sleptons; trilinear couplings among the scalars 
(proportional to the couplings in the superpotential); and gaugino masses for the bino (Mi), 
the winos (M 2 ) and the gluino (M 3 ). Assuming supersymmetry breaking from a hidden 
sector in minimal supergravity (mSUGRA), the SUSY breaking terms are assumed to be 
universal at the scale of grand unification (near the Planck scale) and denoted as mo, Aq 
and M1/2, respectively. Hence the parameters of the corresponding cNMSSM are, apart 
from the gauge and quark/lepton Yukawa couplings, 

mo , Aq , M1/2 , A and k . (2.2) 

It is convenient to fix k from the requirement that the Higgs vevs h u and hd generate the 
correct value of Mz- 

As mentioned in the introduction and discussed in detail in [ZD [22], the remaining 
parameters in ( 12. 2 p are strongly constrained: mo must be small such that the vev s is non- 
zero. A small non-zero value for mo affects essentially only the singlet-like CP-even Higgs 
mass [22] , which is irrelevant for the present study; hence we assume m = in the following. 
In order to avoid a t\ LSP, the singlino (the fermionic component of S) must be lighter such 
that T\ is the NLSP. The singlino relic density can be reduced to an amount compatible 
with WMAP constraints, if its co-annihilation rate with the t\ is large enough, i.e. if the 
corresponding mass difference is sufficiently small. This fixes Aq ~ — \M\/2 [2Tl[22] . Finally 
A must also be quite small, since A induces mixings in the CP-even Higgs sector between 
the doublet- and singlet-like Higgs states: if the singlet-like Higgs state is heavier than the 
(Standard Model like) Higgs state h, the mass of h falls below the LEP bound of ~ 114 GeV 
if A is too large; if the mass of singlet-like Higgs state is below m^ (below 114 GeV as for the 
point P520 in Table [1] below) , its coupling to the Z boson violates again LEP bounds [M] 

for A too large. All in all one finds A < 0.02 (2HE2] ( but A ~ 10 ~ 5 in order stiU to allow 
for singlino-ri co-annihilation). 

A induces also mixings between the singlet-like neutralino and the MSSM-like neutralinos 
(bino, neutral wino and higgsinos). For A < 0.02, these mixings are very small. Hence the 
couplings of the singlino-like LSP Xi to all MSSM-like sparticles (squarks, gluino, sleptons, 
charginos and neutralinos), which are induced by these mixings, are very small as well. 
Accordingly branching ratios of all these sparticles into the singlino-like LSP are negligibly 
small, unless a decay into Xi is the only decay possible. Due to R-parity conservation this 
is the case for the NLSP, the T\. Hence sparticle decay cascades proceed as in the MSSM 
(with a spectrum as in the cNMSSM, but without the singlet-like states) until the T\ NLSP 
is produced. Depending on A, the width of the final decay T\ — > Xi + r can be so small, 
that the T\ decay vertex is visibly displaced [22]. This case (where the displaced vertex 
corresponds to the production of a soft r-lepton) could be another interesting signature for 
the cNMSSM, but subsequently we will not assume that A is so small that this phenomenon 
occurs. 
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Concerning the remaining parameter M1/2, we find that the LEP constraints on the 
Higgs sector require Mi/ 2 > 520 GeV. Then, all bounds on sparticle masses from colliders 
as well as constraints from B-physics are satisfied. For the calculation of the spectrum 
we use the code NMSPEC [23] within NMSSMTools [35l|36], updated including radiative 
corrections to the Higgs sector from [37]. The dark matter relic density is computed with 
the help of micrOMEGAs [9]. Clearly, very large values of Mi/ 2 are generally disfavoured 
by fine-tuning arguments; moreover, smaller values of Mi/ 2 < 1 TeV allow to explain the 
discrepancy of the measured anomalous magnetic moment of the muon with the Standard 
Model [22]. Subsequently we confine ourselves to Mi/ 2 < 1 TeV. In Tabled] we show the 
values for Aq, tan/3 = h u /hd and /i e // (Aq is determined by the correct dark matter relic 
density, whereas tan/3 and fi e ff are obtained as output) as well as the Higgs and sparticle 
spectra for M 1/2 = 520, 600, 800 GeV and 1 TeV for m = and A = 0.001. 

We see that, as announced above, — m x o < 5 GeV and hence r-leptons from the 
decay T\ — > Xi + r are necessarily soft, r-leptons from the decay x 2 ~~ * T\ + r (where x 2 is 
dominantly bino-like) profit at least from m x o — ~ 70 GeV (for P520), or more energy 
from the decays of other sparticles into T\. Note that right-handed sleptons e# and JIr 
decay essentially via the three-body channel as ejj — >■ e + t\ + r. In fact, apart from the 
NMSSM-specific decay T\ — > Xi + T (with a branching ratio of 100 %), the corresponding 
sparticle decay branching ratios can be obtained from the code SUSY-HIT [38] and the 
MSSM with a corresponding spectrum, which is used for the simulations of events below. 

At the LHC, the dominant sparticle production processes are of course squark-gluino 
and squark-(anti-)squark pair productions. Subsequently a typical squark decay cascade 
looks like 

q^q + xl^q + r + n^q + r + r + Xi, (2-3) 
but many more possibilities exist. Their simulation, together with the simulation of Stan- 
dard Model background processes, will be discussed in the next sections. 

3 The signal, backgrounds and cuts 

As for most SUSY models, the production of squarks of the first generation and of gluinos 
will be the dominant sparticle production processes at the LHC. Their total production 
cross sections are obtained at NLO (QCD) from PROSPINO [39], and are also shown in 
Table [TJ The dominant contributions originate from squark + gluino production (~ 37%) 
and squark pair production (~ 25%); less dominant are squark + antisquark production 
(~ 13%) and gluino pair production (~ 5%). The production cross sections of stop and 
sbottom squarks, sleptons, charginos and neutralinos add another ~ 19% to the total 
sparticle production cross sections. 

The dominant background processes for SUSY searches are well-known: top-antitop 
pair production, W+n-jet production, Z+n-jet production, W+Z production and WW+n- 
jet production. Since we will compare the performance of our simulation with the results 
of standard SUSY searches by ATLAS [10] , we assume the same production cross sections 
for these background processes as in [10] (given in Table [2] below). 

Given that gluinos (whose decay generates typically two hard jets) are even somewhat 
heavier than the first generation squarks (generating typically one hard jet), we require at 
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Table 1: Input parameters, tan/3, /i e jj and low-energy spectra for four points of the 
cNMSSM with m = and A = 0.001. In the last line we give the total NLO cross 
sections for the production of all sparticles at the LHC. 

least two hard jets per event only whereas in generic supersymmetric analyses usually four 
hard jets are required. On the other hand, given the large squark and gluino masses, we 
can require quite large transverse momenta of the jets as well as a large missing transverse 
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energy E!f lss . 

For the r-leptons we consider their hadronic decays only. For their transverse momenta 
we require at least 30 GeV, which allows to assume an efficiency of ~ 40% [10] (and a r-fake 
rate of jets of ~ 1 — 2%). On the other hand, since only two among the four r-leptons per 
event are sufficiently energetic and the total signal cross sections are already quite small, 
we require one identified r-lepton only. 

Additional standard cuts are a lower limit on the angle A$ between the hard jets and 
E™ ms , as well as a cut on the transverse mass formed from E™ 1SS and the identified 
r-lepton (in order remove semileptonically decaying W+jets events). Altogether, the list of 
our cuts is given by: 

1. At least two jets, one with p^ > 300 GeV and one with pp > 150 GeV 

2. E^ iss > 300 GeV 

3. At least one r-lepton with p^ > 30 GeV 

4. A$(j i7 E^ s ) > 0.2 for the hard jets 

5. Mj> > 100 GeV, where My is computed from the visible momenta of the hardest 
r-lepton and E™ 1SS . 

Below we will denote this set of cuts as cNMSSM analysis. 

4 Simulation, signal and background rates 

Both the signal and the top quark background were generated by PYTHIA 6.4 [41], which 
was in charge of generation and phase-space decays. PYTHIA performed the parton show- 
ering as well as the matching procedure according to the MLM prescription including initial 
and final state radiation. TAUOLA [4"3TT4"o] was employed for the r-decays. All other back- 
grounds (involving at least one W or Z-boson) were generated with ALPGEN [12]. In order 
to keep the statistics manageable, preselection cuts were applied on the ALPGEN sam- 
ples. Since event generation was performed by leading order generators, the cross sections 
were scaled according to the NLO cross sections in Table [1] for the signal, as in [40] to the 
NLO+NLL calculation for top production, and to NLO (or NNLO level, where available) 
for electroweak boson(s) production. 

For the detector simulation we employed AcerDet [IB]. AcerDet is a fast detector sim- 
ulation which provides a reasonable description of the performance of an LHC detector. 
The events generated by PYTHIA and ALPGEN+PYTHIA were all passed through Ac- 
erDet. AcerDet reconstructs jets, leptons and the missing transverse energy, it also labels 
the origin of the jets, e.g., those coming from a tau lepton. The efficiency and the corre- 
sponding background rejection for a working point of 40% r identification efficiency were 
implemented at reconstruction level. 

One of the issues to be checked is the energy of the reconstructed tau-jets. For this 
initial check we used PYTHIA to produce Z bosons and their subsequent decay to tau 
leptons. The hadron-hadron as well as the lepton-hadron final states were reconstructed 
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and compared to the results of the ATLAS collaboration presented in [3D]. Reasonable 
agreement at the level of a few percent (about 3 GeV) was found. 

A fast detector simulation as AcerDet will not be able to simulate the non-Gaussian 
tails, e.g., in the transverse missing energy distribution. Simulation and reconstruction 
results without full detector simulation and reconstruction are therefore not expected to 
be perfectly reproduced. To get a feeling of how well the background can be modeled with 
AcerDet, two signatures of Ref. [40] have been implemented and analysed in addition to the 
dedicated cNMSSM analysis: the four-jet SUSY search (4j01) and the SUSY search with at 
least one tau in the final state (4jtau). 
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Table 2: The number of simulated background events, cross sections before cuts, and cross 
sections after the 4j01, 4jtau and the cNMSSM analysis. The quoted error is the statistical 
error. 



In Table |2] the result of the background cross sections before and after the cuts is shown 
for the two ATLAS as well as the cNMSSM signatures. For the 4j01 analyses, the AcerDet- 
result for the total background cross section of about 606 fb after cuts is of the right order of 
magnitude compared to 708 fb (NLO) in |40j. Adding the tau identification to the analysis 
in the 4jtau analysis, a total background cross section of 61 fb after cuts obtained here is to 
be compared with 51 fb at NLO in [40]. Thus with and without tau identification, AcerDet 
provides a reasonable estimate of the expected background with respect to the dedicated 
full simulation. 

We note that already with the softer cuts within the 4j01 analyses and notably the cuts 
on A$(jj, i?™ lss ), the remaining QCD background was found to be small in [40]. In our 
case, QCD events could pass the cNMSSM cuts only if a very large value of E™ 1SS and a 
r-lepton would be faked simultaneously. Assuming a jet — > r fake rate up to ~ 2% (for an 
acceptance of 40%), and that the suppression rate of QCD events without missing energy 
for E™ lss ) > 150 GeV is 1% gDJ while we cut at 300 GeV, we should be safe of the QCD 
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background. 

Whereas the two ATLAS analyses are designed to cover a large variety of supersym- 
metric signals, the signatures discussed above are chosen specifically for heavy squarks and 
gluinos as well as r-rich final states as in the cNMSSM. The background cross sections 
for this analysis are shown in the last column of Table [2j The total background, already 
decreased from the 4-jet-0-lepton to the 4-jet-tau analysis by an order of magnitude, is 
reduced by another factor four to ~ 16 fb. 

Typically the overall efficiency for all SUSY processes weighted by the cross section 
varies between 7% and 10%. The cross section for the cNMSSM benchmark points after 
all cuts are shown in Table [3] for the three analyses. The S/B ratio, the ratio S/a/B for an 
integrated luminosity of 1 fb -1 as well as for 30 fb -1 are shown. 
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Table 3: Signal expectation for the NMSSM points at NLO after all cuts for the benchmark 
points. At least 120000 events per point were generated. The error is the statistical error. 
For every point the ratios S/B and S/vB for an integrated luminosity of 1 fb" 1 and 30 fb -1 
are shown. 

Table [3] clarifies the advantage of the cNMSSM cuts with respect to the general analysis: 
the ratio S/VT3 is 7—10 times larger for the cNMSSM cuts even with respect to the standard 
4-jet-tau signature, which originates both from the larger background suppression and the 
larger efficiency on the signal. Correspondingly the cNMSSM cuts allow for a sensitivity, 
for a given luminosity, on a much larger part of the cNMSSM parameter space (for heavier 
squarks/gluinos). The point P800 (with squark/gluino masses of 1.6/1.7 TeV) is hardly 
visible within the standard analysis even for 30 fb -1 , whereas the ratio S/a/B is still ~ 18 for 
the cNMSSM cuts. Only for the point P1000 (with squark/gluino masses of 1.9/2.2 TeV) 
a larger luminosity and/or even harder cuts seem to be required for detection. 
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In Figure [T] we show the spectrum of the effective mass M e g = p£ ts + ^ Pj? + E!f lss 
after all cuts of the cNMSSM analysis are applied, normalised to an integrated luminosity 
of 1 fb _1 . Typically, the spectrum of the effective mass peaks at a value corresponding to 
the masses of the pair produced sparticles [ID]. Here the maxima of M g are shifted to 

'miss 
T 

M e g is harder for the points with heavier squarks/gluinos. 



somewhat larger values due to the cuts on p£ s and E™ 1SS . As expected, the spectrum of 




2000 2500 



3000 3500 4000 
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Figure 1: The effective mass distribution for the SM background and the cNMSSM points 
from Table 1, after all cuts of the cNMSSM analysis are applied, normalised to an integrated 
luminosity of 1 fb _1 . 

For completeness we show in Figure [2] the transverse momentum of the leading tau 
candidate (after the cNMSSM cuts). Modulo the rate, the spectrum of the leading tau 
candidate is slightly harder for the points with heavier squarks/gluinos. 



5 Comparison with the cMSSM in the stau coannihi- 
lation region 

It is well-known that r-rich final states from squark or gluino production would also be 
generated in the so-called stau coannihilation region of the (c)MSSM [29TI33] . Hence the 
question arises, by which signatures this region of the cMSSM can be distinguished from 
the cNMSSM. 
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Figure 2: The transverse momentum distribution of the leading tau after all cuts of the 
cNMSSM analysis are applied, normalised to an integrated luminosity of 1 fb _1 . 



Clearly, the neutralino/r spectrum of the cMSSM is different from the one of the cN- 
MSSM: in the cMSSM (in the stau coannihilation region, which we assume henceforth) the 
lighter T\ has a mass close to the bino-like neutralino LSP x?> whereas the neutralino x\ 
is typically wino-like. At first sight, a squark decay cascade as in Eq. fj2.3j) is also possible 
within the cMSSM, with corresponding replacements of the neutralinos Xi an d X°- How- 
ever, all right-handed squarks (and sleptons) would not couple to the wino-like X21 an d 
prefer to decay directly into the bino-like Xi- These decays do not lead to two r-leptons in 
the cascade. As a consequence, the r-rich cascades hardly occur for right-handed squark 
decays, and are thus less frequent (relative to the total squark production cross section) 
than in the cNMSSM. 

On the other hand, squarks and gluinos can be considerably lighter in the cMSSM than 
in the cNMSSM, since smaller values of Mi/ 2 - together with larger non-zero values of mo 
- are allowed. In particular this is the case, if we look for a point in the cMSSM parameter 
space with similar Xi an d f\ masses as the point P520 of the cNMSSM. It turns out that, 
for the corresponding values of M1/2 and m (we take A = for simplicity), the squarks 
and gluinos are considerably lighter than for the point P520. Hence we denote this point 
as MSSM1 ("1" for "light"). Its parameters and sparticle masses are given in Table H] below. 

As indicated in the last line in Table HI the lighter squarks and gluinos imply considerably 
larger production cross sections for the point MSSM1 compared to P520. As a consequence, 
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Table 4: SUSY breaking parameters, tan/3, fi e g, sparticle spectra and total sparticle cross 
sections for the cMSSM points MSSM1 and MSSMh. 



the number of events passing our cNMSSM analysis above is larger than for P520, in spite 
of the absence of r-leptons in right-handed squark decays. 

There exist also points in the cMSSM parameter space where the squark and gluino 
spectrum resembles the one of P520, implying similar production cross sections. Such points 
correspond to larger values of M1/2 and mo; an example is given by the point MSSMh ("h" 
for "heavy"), whose squark and gluino masses are similar to those of the cNMSSM point 
P520 (see Table gj. 

The signal rates after the cNMSSM cuts for the points MSSM1 and MSSMh are given 
in Table |5j these are considerably larger (as compared to P520) for the point MSSM1, but 
smaller for the point MSSMh in spite of the similar squark/gluino masses and hence the 
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similar total sparticle cross section (see Table H]). The reason was mentioned above: right- 
handed squarks do not decay via r-rich cascades and, hence, right-handed squark decays 
do not contribute to the signal after the cNMSSM analysis. 
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Table 5: Signal expectation for the MSSM points at NLO after all cuts. 120000 events per 
point minimum were generated. The error is the statistical error. 




Figure 3: The effective mass distribution for the points P520, P600, MSSM1 and MSSMh 
after all cuts of the cNMSSM analysis are applied. 

The M e ff spectrum for the MSSM points, together with P520, P600 and the SM back- 
ground, is shown in Fig. El First, the point MSSM1 (with similar \i an d f\ masses as the 
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point P520) has not only a larger signal cross section as compared to P520, but we see that 
its maximum of M e g is visibly shifted towards smaller values. 

Can we distinguish the point MSSMh from any of the cNMSSM points? Due to the 
similar squark/gluino masses as P520, MSSMh has its maximum of M c g in the same region 
as P520, but a significantly smaller signal rate. The signal rate for the cNMSSM point P600 
of about 58 fb is still larger than 41 fb for MSSMh. (The difference is slightly larger than a 
conservative error of 20-30% on the theoretical cross section prediction.) On the other hand 
we see in Fig. [3] that the maximum of M e g- for P600 is shifted towards larger values due to 
the heavier squarks/gluinos: the root mean square of the distributions is about 500 GeV 
and the difference of the average effective mass is about 130 GeV, so that the error on the 
average effective mass for 1 fb -1 is about 70 GeV, i.e., about two times smaller than the 
difference. The average effective mass is affected somewhat by the tails for large effective 
mass. Using a simple fit of the distributions, the peak to peak difference increases slightly 
to about 150 GeV providing for a stronger separation of MSSMh and P600. Any cNMSSM 
point with still heavier squarks/gluinos (such that the signal rate coincides with the one for 
MSSMh) will imply a maximum for still larger values of M e g. Hence, the cNMSSM points 
have either measurably larger signal rates after applying the cNMSSM cuts (if the maxima 
of M e g coincide with a MSSM point), or maxima at measurably larger values of M e g- (if the 
signal rates coincide with a MSSM point). 

Additionally one can compare the cross section ratios after the generic supersymmetric 
4j01 cut, which are 242 fb (MSSMh, where squarks are somewhat lighter than gluinos) 
as compared to 47 fb (P600, where gluinos are somewhat heavier than squarks). Hence, 
given a corresponding signal in the data, a careful comparison of both the signal rates for 
generic and dedicated searches and the maximum of M e g (possibly including in addition 
the transverse momentum of the tau lepton) should allow to distinguish the cNMSSM from 
the MSSM in the stau coannihilation region. 

6 Conclusions and outlook 

In the present paper we have proposed criteria for the search for the fully constrained 
NMSSM at the LHC. In view of the relatively heavy squarks and gluinos in the cNMSSM 
and correspondingly small production cross sections, this task is not quite trivial. On the 
other hand, due to the large number of r-leptons in the final states, signatures involving 
hadronic r decays are relatively efficient. Whereas the soft r-leptons in the final states are 
difficult to use, the requirement of at least one hard r-lepton has a relatively large signal 
acceptance. 

Combining this requirement with relatively hard cuts on the transverse momenta of two 
jets and E™ ss as specified at the end of section 3, the signal to background ratio is signifi- 
cantly improved with respect to the more standard 4j01 or 4jtau analyses. This result was 
obtained after simulations including detector effects and a r acceptance, which we com- 
pared with and checked against the analysis of SUSY signals by the ATLAS group. Hence, 
already an integrated luminosity of 1 fb -1 (at 14 TeV cm. energy) becomes sensitive to 
part of the parameter space of the cNMSSM whereas, trivially, more luminosity is required 
in case of heavier squarks and gluinos. In any case we believe that the cuts proposed here 
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are the most sensitive ones to the parameter space of the cNMSSM. In addition we have 
discussed in how far a refined analysis employing both the signal rate and the maximum of 
M e fj allows to distinguish the cNMSSM from the MSSM in the stau coannihilation region. 

In the near future the LHC is on track to accumulate an integrated luminosity of 1 fb _1 
at 7 TeV cm. energy at the end of 2011. We have estimated the number of signal events 
for the point P520, if we lower the cNMSSM cuts correspondingly: for two jets we require 
Pt > 50 and 20 GeV, respectively, and E™ 1SS > 200 GeV. We obtain about 5 signal events 
passing these cuts. 

For our analysis of signals of the cNMSSM at 14 TeV cm. energy we have left aside 
the presence of two soft r leptons per event which represent, in principle, a spectacular 
signature for this class of models. In some regions of the parameter space of the cNMSSM 
- for very small values of A and/or a small T\-x\ mass difference - the life time of T\ can be 
very small leading to displaced vertices of the decay r x — > r + \\ into these soft r-leptons. 
Using dedicated track-based algorithms, the search for these soft r-leptons originating from 
displaced vertices is perhaps not completely hopeless. 

In the framework of the general NMSSM, a singlino-like LSP can be accompagnied by a 
NLSP different from t\. The signatures of these scenarios would be very different from the 
ones discussed here (depending on the nature of the NLSP), and should also be investigated 
in the future. 
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